Twinlike Models in Scalar Field Theories 
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This work deals with the presence of defect structures in models described by real scalar field in a 
diversity of scenarios. The defect structures which we consider are static solutions of the equations 
of motion which depend on a single spatial dimension. We search for different models, which support 
the same defect solution, with the very same energy density. We work in flat spacetime, where we 
introduce and investigate a new class of models. We also work in curved spacetime, within the 
braneworld context, with a single extra dimension of infinite extent, and there we show how the 
brane is formed from the static field configuration. 

PACS numbers: 11.27.+d, ll.10.Kk 



I. INTRODUCTION 

Topological structures appear in several areas of non- 
linear science. In high energy physics, the most known 
topological structures are kinks, vortices and monopoles. 
Kinks are the simplest structures, and they appear in 
one spatial dimension, in the presence of a single real 
scalar field. Vortices and monopoles are more involved, 
and they usually appear in two and three spatial dimen- 
sions, in the presence of complex scalar fields coupled to 
Abelian and non Abelian gauge fields, respectively [H, 0- 

In the present work we focus our attention on kinks 
in relativistic systems described by real scalar field. The 
subject has been investigated in a diversity of contexts, 
but here we bring novelties related to the issue recently 
studied by Andrews, Lewandowski, Trodden and Wes- 
ley (ALTW) in Ref. Q, concerning the presence of twin 
models, that is, of distinct models which support the 
very same topological structure. To make the investi- 
gation more general, however, we will consider models 
engendering topological or kinklike and nontopological 
or lumplike structures. 

Relativistic models described by scalar fields have been 
studied in a diversity of scenarios. Here we quote some 
works dealing with junctions of kinklike structures 041] i 
the deformation procedure which lead us to new models 
and the respective kinklike solutions [1, [l(| , the presence 
of vortons [ll[ and other new models and solutions fl2i — 
[13]. Also, from the experimental point of view, new and 
interesting kinklike solutions have been found recently in 
systems constrained to work within small regions, at the 
nanometer scale [H, [l9| . 

Motivated by cosmology, another direction of study 
concerns models where kinematic modifications of the 
scalar field are introduced, attempting to contribute to 
explain the present accelerated expansion of the Universe 
[2(j-[23[ . Recent investigations of scalar fields with modi- 
fied kinematics have been presented in f24T - l30l ] with dis- 
tinct motivations, in particular to search for the presence 



of compacton, that is, of topological structures which 
live in a compact region, and for applications in the 
braneworld context, within the five dimensional space- 
time engend ering a single extra dimension of infinite ex- 
tent |31H33| |. 

In the study of models with modified kinematics, in 
3] the authors noted that it is possible for k-defects to 
masquerade as canonical scalar field solutions. That is, 
given a standard scalar field model described by the field 
4>, with canonical kinetic term and potential V((f>) in the 
usual form, the topological defect profile (p(x) and cor- 
responding energy density p(x) can also appear in a k- 
defect theory. We call the modified or k-defect theory 
introduced in H, the ALTW model. And we recall that, 
despite having identical defect solutions, the two theories 
are not reparametrizations of each other, since the fluctu- 
ation spectra about the kink are different. In the current 
work, we call the two distinct models twinlike theories, 
because they support the very same defect structures, 
engendering the very same energy density. 

In this work we further investigate this issue, consid- 
ering as twinlike theories the standard model and the 
ALTW model introduced in [3]. We deal with the de- 
fect solutions, and we compare the two models in the 
case where the static solutions are more general. We 
do this studying solutions of the equations of motion in 
the (4>, 4> r ) plane, where prime stands for spatial deriva- 
tive of the field, 0' = d<fi/dx, with <f>(x) describing static 
solution. Also, we extend the results to the case of non- 
topological or bell-shape solutions, which are of direct 
interest to tachyon condensation in string theory [34- 37], 
and to cosmology [3o |. 

We take advantage of the formalism introduced in 
Ref. [29j], from which it is possible to get to first-order 
differential equation that solves the equation of motion 
for defect solutions of a general scalar field model. Wc 
use this to verify if the issue is specific of the two models 
introduced in Q. In the process, we further ask for other 
models, with the answer offering a new class of twinlike 
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models. We also study the issue in curved spacetime, cou- 
pling the scalar field with gravity within the braneworld 
context, with a single extra dimension of spatial nature 
and infinite extent. 

To ease understanding of the subject, we organize the 
work as follows: In the next Sec. [IT] we explore the stan- 
dard and the ALTW models, and we comment on some 
specific properties the two models engender. In Sec. Hill 
we investigate two distinct ways to distinguish twinlike 
models. In Sec. IIVI we generalize the investigation to an- 
other class of twinlike models, and in Sec. [V] we deal with 
twinlike models in the braneworld context. We end the 
work in Sec. I VII where we include comments and conclu- 
sions. 



II. THE MODELS 

Let us now focus attention on the standard and mod- 
ified models. We first consider the standard model, and 
then we deal with the modified model. We also comment 
on some specific properties the two models engender. 

A. Standard model 

The standard model is described by the Lagrange den- 
sity 

= x-v{4>), (i) 

where 

x = -dpm, (2) 

and V{4>) is the potential. The equation of motion has 
the form 

d^(j> + V+ = 0, (3) 

where stands for the derivative of the potential with 
respect to the field 0, that is, V$ — dV/dcf). The energy- 
momentum tensor has the form 

= d»(f>d^-7f£ s . (4) 

Here we work in (1, 1) space-time dimensions, with the 
metric obeying diag^" = (1,-1). The equation of mo- 
tion for static solution becomes 

<t>" = V<f, (5) 

where 0' = d<p/dx, etc. And the energy-momentum ten- 
sor reduces to 

T 00 =p(x) = ±<t>' 2 + V(<P), (6) 



which is the energy density, and 

T 11 — ^0' 2 ~ V{4>)- (7) 

As an example, we consider the potential with quartic 
self-interaction, in the form 

^) = i(W 2 ) 2 , (8) 

taking the field and coordinates as dimensionless quan- 
tities. This model is well known, and it has the minima 
<f>± = ±1. Also, it supports the defect solutions 

4>±{x) = ± tanh(x — x), (9) 

where x stands for the center of the defect, which we 
take to be at the origin in the x-axis, that is, x = 0. The 
two minima 0± define a topological sector, and the above 
solutions are topological or kinklike defects. 
We can consider another model, with potential 

V sG {4>) = \ cos 2 (0). (10) 

This is the sine-Gordon model, and typical solutions are 

<t>{x) = ± arcsin(tanh(a:)) + kn, k = 0, ±1, ±2, . . . (II) 

We can use another potential, with cubic self- 
interaction, in the form 

V3W>H^ 2 (W). (12) 

This model has a single, local minimum, at <fi = 0, and it 
supports the defect solution 

4>{x) = sech 2 (a;). (13) 

This solution describes a nontopological or lumplike de- 
fect. Note that (Ti"2"j) has a partner potential, with the 
factor (1 — 0) changed to (1 + 0); the partner model sup- 
ports the non-topological defect <j>{x) = — sech 2 (x). 
Another model is described by the potential 

VmW = ^ 2 (1 - <P 2 )- (14) 

It is an inverted <fi 4 model, and it has a local minimum at 
= 0; it supports the nontopological or lumplike defects 

<j>{x) = ±sech(ir). (15) 

B. The ALTW model 

Let us now turn our attention to the modified model 
which was considered in Ref. [3J. It is described by the 
Lagrange density 

£,, = M*-M*(l+W>)0Tf. {16) 
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We first note that in the limit where the mass scale M 2 
is very large, wc get from the modified Lagrange density, 
up to the order 1/M 2 : 



x - u(4>) + 



i f i 



M 2 \ 2 



-X 2 + U((f>)X 



(17) 



Thus, we see that the limit 1/M 2 — > 0, leads us back to 
the standard model ([T]), with the identification U(4>) — 
V (</>), so we call U (4>) the potential of the modified model. 

In general, for the modified model (fT6"| . the equation 
of motion is given by 



1 



u{4>) 
~TP~ 



1 



IX 
TP 



2X „ 

m 2 = 



(18) 



or more explicitly 

1 d^cj>d v <i 



7 



M 2 l-^ 



d^d u 4> = - 



M 2 i 



u(<t>) 

M 2 



(19) 



Here the energy-momentum tensor has the form 
(1 + ^)5^0 



rpfll> _ 



2X 

TP 



M 2 rj^ 1- 



§i< 2 »» 



We focus our attention on static configurations, with 
= cj>(x); in this case, the equation of motion becomes 



p{x) 



-M + M 



U{4>) 
TP" 



-M 



1 - 

1 TP 

c 



A'2 



M U-Ip ) (1 + T72 )-(24b) 



Af 2 



In Ref. [25J one shows that a necessary condition to have 
stable kinklike solutions is to set C = 0. In this case, the 
equation ((23)) becomes 



4> a = 2U{4>) + 



U 2 (4>) 

M 2 : 



and the energy density gets to the form 

u 2 (<t>) x 



p{x) = 2t/(0) 



M 2 



(25) 



(26) 



C. Some specific properties 

Let us now follow |3j and consider the standard model, 
with the specific potential 



(27) 



where U (0) is the potential for the modified model. We 
note in the standard model with the above potential, the 
equation of motion is given by 



rrdnW + 1 + 



M 2 



u+ = o, 



(28) 



1 



1 + M 2 1 

1 + M 2 1 



U(S) 

~TP~ 



(21) 



We integrate this equation once, to get the first-order 
equation 



1 - 



1 



U(6) 

~TP~ 



1 



C 

M 2 ' 



(22) 



where C is an integration constant. This equation can 
be rewritten as 



A'2 



-M + M 2 



1 



EM 

A/ 2 



1 



C 
M 2 



(23) 



Since the energy-momentum tensor is conserved, for 
static field we get that T 11 has to be a constant, which 
leads us to Eq. (|2"2")l . On the other hand, the energy 
density T 00 = p(x) has the form 



'"■'■' = - 1/2 + 1/2 ( 1 + i^)V 1+ w> (24a) 



and the energy-momentum tensor has the form 



1 U 2 ^) 



= d»4>d v <\> - rf v X + rf v U{4>) + - 



2 M 2 



In the case of a static solution, we get 

U(<f>) 



and 



= 1 



M 2 



P{x) = ^ 2 + U{4>) + -^^. (31) 



Us, 



l[/ 2 («/>) 



However, from (|2T)1) we can write 



(29) 



(30) 



t 11 = ^4>' 2 — U{4>) — ( 32 ) 

Since the energy-momentum is conserved, for static so- 
lution we can write T 11 — C, and this leads us to the 
equation 



(j) 12 = 2C + 2U(4>) + 



u 2 (<t>) 

M 2 ' 



(33) 
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and for C = one gets to 



cf>' 2 = 2U((f>) + 



u 2 (<t>) 

M 2 ' 



(34) 



We see that solutions of the above first-order equation 
solve the equation of motion (|30j) . Thus, we reproduce 
the result of Ref. [3j, showing that the modified model 
(flo) and the model (|27p give the same defect structure. 
Moreover, using (|3"4l in (j3"Tj) we reproduce (j2"6"1) . and so 
the defect has the very same energy density, as also stated 
in Ref. 0. 

The above result requires that we make T 11 = 0. In 
the more general case, we can use T 11 = C, and this 
makes the two models different. Therefore, a nice way to 
distinguish between the two theories is simply to inves- 
tigate the solutions taking T 11 = C, as a constant but 
nonvanishing quantity. We explore this fact in the next 
Sec. ED 

We now take advantage of the result firstly presented 
in Ref. [29j , from which one can write a first-order frame- 
work for defect solutions of a general scalar field model. 
The approach concerns introducing the function W (((>), 
such that 



(35) 



for some generic Lagrange density C. In the case of the 
standard Lagrange density, we have L s x = 1, and this 
leads to the equation 



4>' = W h 



It is a total derivative, so the energy follows easily, in the 
form E = |AW| = \W{<j>{po) - W (</>(- oc))\. 

A nice information follows when one expands (|39p for 
small 1/M 2 ; up to the lowest order we get 



u(4>) 



l. 



(41) 



(36) 



On the other hand, from (|2"T)l. the limit 1/Af 2 -> leads 
us to V(4>) = £/(</>). Thus, we have the same W = W{<j>) 
controlling both the standard and modified models, al- 
though V{4>) and U{<f>) are different functions of the real 
scalar field. 



III. DISTINGUISHING TWIN THEORIES 

Let us now focus on the issue of distinguishing the stan- 
dard and modified models. In Q, the authors showed 
that the study of linear stability of the defect structure 
is able to distinguish the two models. Since linear stabil- 
ity is important for introducing quantum corrections at 
the one-loop level, we can say that the quantum effects 
distinguish the two models. More than that, we have 
already shown that there is another way to distinguish 
the two models, if one deals with a constant but non 
vanishing T 11 ^ 0. Thus, below we investigate the two 
possibilities, of using a constant but nozero T 11 , and of 
studying linear stability. However, since linear stabilit; 
was already investigated in Q, here we depart from 
and follow an alternative way to investigate linear stabil- 
ity. 



such that 



v{4>) = 2^1- 



(37) 



This is the first-order framework one usually gets for 
scalar field with standard dynamics. 

However, the very same framework can be used to 
study the modified model. In this case, taking L m x4>' — 
Wa, leads to the equation 



A'2 



urn 

M 2 



W 2 
M 2 



To make this equation compatible with 
we have to write 



(38) 



for C = 0, 



U(4>) 



-M 2 + M 2 Wl + ^. (39) 



An important point in the above procedure is that the 
energy density, which can be written as p(x) — 4> , leads 
us to the nice result 



p(x) = <A 



dW 
dx 



(40) 



A. Nonzero T 

Let us firstly consider the case of a non vanishing T . 
In this case, the equation of motion for the standard and 
modified models are, respectively, 



= Wi 



= -M l 



2C, 



M 



1 



M 2 



(42a) 
(42b) 



We take as an explicit example, the model described by 

W^ = l-<j) 2 . (43) 

Note that for the standard model, the potential V{4>) 
leads to the (/> 4 model, with spontaneous symmetry break- 
ing; see Eq. d8]). With this choice, the first order equa- 
tions (|42a|) and (|42bl) become 



= (1- 
= -1 



r y + 2c, 

1 "T M 2 



(44a) 
(44b) 



for the standard and modified models, respectively. We 
plot in Fig. 1 the solutions in the ((f), <j)') plane, for C > 0, 
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FIG. 1: Plots of the solutions <j>(x) of the standard (left panel) 
and ALTW (right panel) models, for l/M 2 = 1, 0.5, 0.1, from 
the upper panel to the lower one, respectively. The solid 
(black) lines stand for T 11 = 0, the dashed (red) lines for 
T 11 > 0, and the dot-dashed (blue) lines for T 11 < 0. 
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FIG. 2: Plots of the solutions of the sine-Gordon V s g{4>) (left 
panel) and modified (right panel) models, for l/M 2 = 0.5 and 
0.1; conventions as in Fig. 1. 



motion become 



= cos 2 (0)+2C, 



2 (Q) 



-1 + 



M 2 



(1 



■Y 2 



(46) 
(47) 



These equations are investigated in Fig. 2, and there we 
show how the solutions behave for two distinct values 
of l/M 2 . We note that, in terms of l/M 2 , the models 
behave in a way similar to the previous case. 

We can also investigate models that support nontopo- 
logical configurations. In this case we follow [39| and we 
choose 



C = 0, and C < 0, with dashed (red), solid (black) and 
dotted (blue) lines, respectively. Here we note that only 
the black curves, which correspond to kink and antikink, 
are identical in both plots. We note that the two phase 
diagrams tend to become the same as l/M 2 becomes 
greater and greater. The solid (black) curves correspond 
to the topological kink and antikink configurations 



(x) = ± tanh(x) 



(45) 



The standard model is known as the 
get to the equations 



(48) 

b 3 model. Here we 



<^ 2 (1 



-1 



- <t>) + 2C, 
1 <t> 2 (i-4>) 

(l-l^) 2 



(49) 
(50) 



We can investigate another model, with — cos(<fi). 
In the case of the standard Lagrange density, this corre- 
sponds to the sine-Gordon model. Here the equations of 



These equations are investigated in Fig. 3, and there we 
see how the two models behave for a single value of M 2 . 
In terms of l/M 2 , the two models also behave in a way 
similar to the two previous cases. 
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FIG. 3: Plots of the solutions of the standard V^{(j)) (left 
panel) and modified (right panel) models, for 1/M 2 = 0.5; 
conventions as in Fig. 1. 



B. Linear Stability 

Another way to distinguish the twin models is investi- 
gating linear stability, as done in Ref. [|[ . Here, however, 
we follow another route. We introduce general fluctua- 
tions for the scalar field, 4>{x,t) — <j)(x) + r){x) cos(w t), 
where <fi(x) represents the static solution. For a general 
Lagrange density, we use the fluctuations in the equation 
of motion to obtain up to first order in 77, 

- [(2C XX X + C x ) ?/]' = + 0*W)' + ^ 2 Cx)v- 

We first suppose that the solution is constant and rep- 
resents a minimum v of the potential, that is, <j> = v. It 
obeys W^v) — 0, and so we get 



n 



wi 



(52) 



which is valid for both the standard and modified models. 
Therefore, linear stability at the minimum state v does 
not distinguish the two models. 

On the other hand, for static solution <j) — (f>(x), it 
is possible to distinguish the two twin theories. In the 
standard theory, we have the following equation for the 
fluctuation 



T)" + V q (x) T) = UJ 2 7], 



where 



V q {x) = W 2 



(53) 



(54) 



is the potential of the Schrodingerlike equation (1531) . 
However, in the modified ALTW model, the equation 
which controls stability is much more complicated. In- 
deed, it has the general form 



(55) 



L0 2 rj. 




We see that this last equation is different from the previ- 
ous one. However, instead of following the investigation 
done in Q, here we introduce new variables, following 
the procedure introduced in Ref. [29|. We change x — > z 
and r\ — > u using the expressions 



dx 




dz, 



which allow writing the Schorodingerlike equation 



U zz + U q (z) U = L0 2 U, 



where 



U q (z) 



(56a) 
(56b) 

(57) 
(58) 



It is hard to write an analytical expression x = x(z), for 
M generic. However, we can consider the case where M 2 
is very large, in order to obtain analytical results in a 
perturbative way. In this case we can write, up to the 
order 1/M 2 , 



M 2 ' 



(59) 



with the choice W(<fi(0)) — 0. Moreover, we note that 
a generic function of the scalar field <fi, F(cj)(x)), can be 
expressed in the form, up to the order 1/M 2 , 



F(cf>(x)) = F{4>(z)) 



{F^W)(z) 
2M 2 



(60) 



and so the potential can be written in the form 



U q {z) = U°(z) 



aw 2M 2 



(61) 



where U*{z) stands for the the potential of the standard 



model. In the case of W$ = 1 



', we get 



U q (z) = 4-6sech 2 (z) 

+ — [4 sech 2 (z)+14 sech 4 (z)-21 sech 6 (z)]. (62) 



In the general case, we have to rely on numerical in- 
tegration to obtain x — x{z). In Fig. [4] we depict the 
behavior of x = x(z) for some values of 1/M 2 , taking 
— 1 — 4> 2 in ([56^) . Note that in the region out of the 
core of the defect, <j) approaches the minima and is 
small, and we can approximate z w x + c where c is a 
constant, and u ~ rj. 



FIG. 4: Plots of x = x(z) for solid, dashed and dot-dashed 
lines, for 1/M 2 = 0.1, 1, 10, respectively. 



FIG. 5: Plots of the potential U q (z) for solid, dashed and 
dot-dashed lines, for 1/M 2 = 0.1, 1, 10, respectively. 



In the general case, the full expression of the potential 
is given by 



U q (z) = U s q {x{z))- 



M 7 



M 2 



2\ 3 



(x(z)) 



3 . 5W? 



"" 1 AM 2 \M 2 



(x(z)).(63) 



7 14 21 28 35 42 49 56 63 70 77 84 91 98 105 



We use this expression to depict the potential in Fig. [5] 
in the case of W$ — 1 — <fr 2 , for some values of 1/M 2 . In 
the figure, we see that the potential starts as the modi- 
fied Posch- Teller potential, and approximates the case of 
a box, for increasing 1/M 2 . We note that although the 
depth of the box does not depend on M, its width in- 
creases as 1/M 2 increases, adding more and more bound 
states into the box. Similar result was also obtained in 
Ref. [3|. In the general case, we can numerically count the 
number of bound states as 1/M 2 increases. The result 
is depicted in Fig. [51 where the dots indicate the value 
of 1/M 2 for which the number of bound states jumps to 
the next one. 



IV. OTHER TWINLIKE MODELS 

Let us now focus on the issue of obtaining other twin- 
like models. Toward this goal, we consider a more general 
Lagrange density, given by l26jj 



on — 1 



V(cf>). 



(64) 



FIG. 6: Number of bound states supported by the potential 
U q (z) as a function of 1/M 2 . 



We note that the case n = 1 leads us back to the standard 
model. For other values of n, we are ch ang ing kinematics, 
and this is of interest in cosmology j20H23j and also in the 
context of compactons, that is, of defect structures which 
lives in a compact region within the real line [24|-[26j . We 
use the more general model to write 



x 



\X 



\n-l 



(65) 



We now follow the procedure introduced in [29| and for 
static solution we get 



V(0 



2n- 1 



IF?- 1 . 

2n * 



with the energy density given by 
p{x) = wf* 



(66a) 
(66b) 

(67) 
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We suppose that the twin theory has the following form 

£ = M 2 + M 2 F(X)G(cf>), (68) 

where M is the mass parameter, as before, and F(X) and 
G((/>) are functions to be determined. For static solutions, 
we can write 



M 2 F X G{4>W = 



1 



2XF X - F 



(69a) 
(69b) 



Using both ((66)) and (|69|) . and considering that the solu- 
tion <j)(x) and W are the same for both theories, we get 
to the following differential equation 

F x (M 2 - 2 n X\X\"- 1 ) = -F\2X\ n -\ (70) 

which can be integrated to give 

F(X) = (M 2 - 2 n X\X\ n - 1 )^ . (71) 
Therefore, the twinlike theory is given by 



£ = M + M—G(<j>) 1-2" 



.X\X 



n— 1\ 2 ? i 



M 2 



(72) 



The function G(4>) in (I69b|) can be written as 



(7 ( , )) = -l/-|l+^^| . (78. 



and the Lagrange density gets to the form 



£ = M 2 -M z 1 



2n 
2n~ 



1 M 2 J \ 



x\x 



n— 1 \ 27i 



M 2 



(74) 

We can use ([59"]) to rewrite the above result in terms of 
U(4>). It has the form 



£ = M — M 



1 - 



x 1-2" 



2n - 
M 2 " 



(75) 



Here we note that the above expression exactly repro- 
duces the ALTW model described by (fTo) , in the limit 
n — > 1. We also note that the energy density is also 
given by (jBTf , so it behaves as in the standard model, as 
expected. 

We can investigate the case in which 1/M 2 is very 
small. Here we get, up to the order 1/M 2 , 



C 



-X\X\ 



1 2n-l 

~AP 8n 2 



V 



2 n X\X\ 



2)i 



2n- 1 



V{!f>)) .(76) 



Thus, the twinlike model ([74]) reproduces the model ([rj4]) 
in the limit 1/M 2 — > 0. This is similar to the previous 
case, so we can say that the two models (|64| and ((74|) are 
twinlike models, and they are generalizations of the pre- 
vious models, described by ([TJ and (fTBj) . It is interesting 
to remark that the above result was obtained within the 
first-order framework put forward in Ref. [29| . 



V. THE BRANEWORLD SCENARIO 

The defect structure which appears in scalar field the- 
ory can be used to mimic thick brane in the braneworld 
scenario with warped geometry and a single extra dimen- 
sion of infinite extent [3lT-l33j|. In this sense, we can ask 
how the twinlike models studied before can be used to 
mimic similar scenario within the braneworld context, in 
the form of a twinlike configuration. 

To investigate this issue, we follow [3(| and we consider 
a five-dimensional action describing gravity coupled to 
the scalar field in the form 



S = J d 5 x^g \--R + £faX) 
where we are using AirG = 1 and 



X = iv A/ ^V 



M , 



(77) 



(78) 



with M, N = 0, 1, 2, 3, 4, running on the five-dimensional 
spacetime. The equation of motion which we obtain for 
the scalar field in this new scenario is given by 



(79) 



(80) 



G /iU V A y B (p + 2XCx4 
where G has the form 

G AB =£ x g AB + CxxV A W B ^ 
The energy-momentum tensor becomes 

Tab = V^V B ^£ X - g AB £ (81) 

The line element of the five-dimensional spacetime can 
be written as ds 2 — e 2A rj^ v dx >J 'dx L ' — dy 2 , where A is 
used to describe the warp factor. We suppose that both 
A and <f> are static, such that they only depend on the 
extra dimension y, that is, A = A(y) and <f> — 4>{y). In 
this case, the equation of motion ((79")) reduces to 

(C x + 2XC X x) 4>" - (2X£ X ^ - £ ) = -4£ X 4>'A', 

(82) 

where prime now denotes derivative with respect to the 
extra dimension. Also, we are using that £x — d£/dX 
and £<p = d£/d(j>, etc. 

We now focus attention on Einstein equations, which 
lead us to the equations 



A" 
A' 2 



-X£ 



X: 



(£ - 2X£ X ) . 



(83a) 
(83b) 
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where X = —<j) l2 /2 for static configuration, as before. If 
we use (JH2), we can write 

(C - 2XC X )' = -4(f>' 2 A' £ x - (84) 
To get to the first-order framework, we suppose that 



A' 



(85) 



and therefore we get the same equation of motion of the 
standard model. Moreover, we can use (|93l) to write 



2U 



El 
M 2 ' 



(95) 



Therefore, the Lagrange density of the twin brane 
model has the following form 



In this case, the equations ()83a[) and (|83b|) leads us to, 
respectively, 



C-2XC X = \w 2 . 



(86a) 
(86b) 



In the case of a standard scalar field, the Lagrange 
density has the form (JTJ, and so the equation of motion 
becomes 



+ 4A'<j>' = 



(87) 



After substituting this result in (|86a[) and (|86L>[) . we get 
to the following first-order equation 



and the potential 



v=\wl-\w\ 



(89) 



Note that the equations ((55)) and ([85} solve equation ([57)1 
for the potential (j89|) . Note also that the energy density 
is given by 



r 



00 



^ ' I l w i _ lw 2 



(90) 



This is the standard situation, but now we have to find 
the twin model. Toward this goal, let us consider a scalar 
field theory governed by the following Lagrange density 



i+ w)r-w + ™> (91) 



C = M 2 - M 2 
where f(4>) is to be determined. We use the equations 



86ap and (|86bp to write, respectively 



1 + M 2 



2X 
M 2 



1 + 



M 2 



2X 
M 7 



(92) 



(93) 



where we have used f(<j>) = W 2 /3, such that, if we sub- 
stitute JMJ) into (JSU) we get to 



(94) 



C = M — M 



M 2 J 



2X 1. 
M 2 



-W 2 



Moreover, the energy density is 



T°° = e 2A [ Z W $ - -W< ). 



(96) 



(97) 



which exactly reproduces the previous expression (|90l) . 
Thus, the two models have the same solution, with the 
very same energy density. This is a feature of twinlike 
models, as expected. In this case, the potential can be 
written as 



U(4>) 




We use this to write 



U{4>) = -M* + M\ 1 1 + 2 



v{4>) 

M 2 



2 W 2 (4>) 
3 M 2 ' 



(99) 



which relates the potentials of the two theories. Here we 
note that, in the limit of very small 1/M 2 , we can use the 
above expression ((99]) to verify that the twinlike model 
(|9U1) lead us back to the standard model, as expected. 



VI. COMMENTS AND CONCLUSIONS 

In this work we extended the idea put forward in 0], 
concerning the presence of two distinct scalar field theo- 
ries which support the same defect structure. We called 
the two field theory models twinlike models, since they 
engender the same defect structure, with the very same 
energy density. 

We have shown that the two models can be distin- 
guished classically, using the fact that T 11 is constant, 
but it does not need to be zero. For nonzero T , the 
two equations of motion give distinct static solutions. 
Distinct results also appear when one investigates lin- 
ear stability of the defect structure in the two models. 
This was firstly presented in [3j, but here we have inves- 
tigated stability in the standard way, searching for the 
Schrodingerlike equation which usually appears in the in- 
vestigation of linear stability. Since the ALTW model is 
described by a modified Lagrange density, of the k-defect 
type, it does not lead us directly to a Schrodinger equa- 
tion. However, we have found the correct way to obtain 
the Schrodinger equation, and this allowed us to search 
for the number of bound states connected to the defect. 
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We have extended the idea of obtaining twinlike models 
to a new class of field theory models. This is interesting 
since it shows that the presence of twinlike models is not 
specific of the pair of models presented in Q. In this 
sense, the presence of twinlike models seems to deserve 
further investigation, since it may be a generic situation 
in Field Theory. 

Another line of investigation concerns the presence of 
twinlike models within the braneworld context [3l|. As 
one knows, scalar fields can be used to generate thick 
branes [32l |33[ , and so we have extended the idea of twin- 
like models to the braneworld scenario. The investigation 
led us to the twinlike braneworld concept, in which a 
brane can be generated from two distinct models of real 
scalar fields, coupled to gravity in (4, 1) dimensions, with 
a single extra dimension of spacelike nature and infinite 
extent. 

The existence of twinlike models poses an interest- 
ing issue, concerning the coupling of scalar fields with 



fermions, which is usually of the Yukawa type. In this 
case, one knows that fermions may fractionalize in the 
presence of topological defect , like the ones we have 
studied in this work. In this sense, although one can con- 
struct two distinct scalar field models to couple to the 
fcrmion field, if the two scalar field models are twinlike 
models, the defect structure they may engender will in- 
duce the very same quantum corrections to the fermions, 
at least at the lowest quantum level. Thus, in the pres- 
ence of fermions, it seems that we need to go beyond the 
first quantum corrections, in order for the fermions to 
distinguish the two models. Another issue refers to the 
coupling with gauge fields, to study the presence of twin- 
like models for vortices and monopoles. These and other 
related issues are presently under consideration, and we 
hope to report on them in the near future. 

We would like to thank CAPES and CNPq, Brazil, 
and FCT project CERN/FP/116358/2010, Portugal, for 
partial financial support. 



A. Vilenkim and E.P.S. Shellard, Cosmic strings and 

other topological defects (Cambridge, Cambridge/UK, [20 

1994). 

N. Manton and P. Sutcliffe, Topological solitons (Cam- [21 
bridge, Cambridge/UK, 2004). 

M. Andrews, M. Lewandowski, M. Trodden, and D. Wes- [22 
ley, Phys. Rev. D 82, 105006 (2010). 

G. W. Gibbons and P. K. Townsend, Phys. Rev. Lett. [23 
83, 1727 (1999). [24 

H. Oda, K. Ito, M. Naganuma, and N. Sakai, Phys. Lett. [25 
B 471, 140 (1999). 
S.M. Carroll, S. Hellerman, and M. Trodden, Phys. Rev. [26 
D 61, 065001 (2000). 

D. Bazeia and F.A. Brito, Phys. Rev. Lett. 84, 1094 [27 
(2000). [28 
D. Bazeia and F.A. Brito, Phys. Rev. D 62, 101701 
(2000). [29 
D. Bazeia, L. Losano, and J.M.C. Malbouisson, Phys. 
Rev. D 66, 101701 (R) (2002). [30 
V.I. Alonso, D. Bazeia, M.A. Gonzalez Leon, L. Losano, 
and J. Mateos Guilarte, Nucl. Phys. 810, 427 (2009). [31 
Y. Lemperiere and E.P.S. Shellard, Phys. Rev. Lett. 91, 
141601 (2003). [32 
D. Bazeia, J. Menezes, and R. Menezes, Phys. Rev. Lett. 
91, 241601 (2003). [33 
M. Gleiser and R.C. Howell, Phys. Rev. E 68, 065203(R) 
(2003). [34 
M. Gleiser, Phys. Lett. B 600, 126 (2004). [35 
G. Basar and G.V. Dunne, Phys. Rev. Lett. 100, 200404 [36 
(2008). [37 
A. Alonso-Izquierdo, M. A. Gonzalez Leon, and J. Ma- [38 
teos Guilarte, Phys. Rev. Lett. 101, 131602 (2008). 
S. Dutta, D.A. Steer, and T. Vachaspati, Phys. Rev. Lett. [39 
101, 121601 (2008). 

P.O. Jubert, R. Allenspach, and A. Bischof, Phys. Rev. [40 
B 69, 220410(R) (2004). 

A. Vanhaverbeke, A. Bischof, and R. Allenspach, Phys. 



Rev. Lett. 101, 107202 (2008). 

C. Armendariz-Picon, T. Damour, and V. F. Mukhanov, 
Phys. Lett. B 458, 209 (1999). 

J. Garriga and V. F. Mukhanov, Phys. Lett. B 458, 219 
(1999). 

T. Chiba, T. Okabe, and M. Yamaguchi, Phys. Rev. D 
62, 023511 (2000). 

R. J. Scherrer, Phys. Rev. Lett. 93, 011301 (2004). 
E. Babichev, Phys. Rev. D 74, 085004 (2006). 

D. Bazeia, L. Losano, R. Menezes and J. C. R. E. 
Oliveira, Eur. Phys. J. C 51, 953 (2007). 

C. Adam, J. Sanchez-Guillen and A. Wereszczynski, J. 
Phys A 40, 13625 (2007); 42, 089801 (E) (2009). 
M. Olechowski, Phys. Rev. D 78, 084036 (2008). 

C. Adam, N. Grandi, P. Klimas, J. Sanchez-Guillen and 
A. Wereszczynski, J. Phys. A 41, 212004 (2008). 

D. Bazeia, L. Losano and R. Menezes, Phys. Lett. B 668, 
246 (2008). 

D. Bazeia, A. R. Gomes, L. Losano and R. Menezes, 
Phys. Lett. B 671, 402 (2009). 

L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 
(1999). 

W.D. Golberger and M.B. Wise, Phys Rev. Lett. 83, 4922 
(1999). 

O. DeWolfe, D.Z. Friedmam, S.S. Gubser and A. Karch, 
Phys. Rev. D 62, 046008 (2000). 

A. Sen, JHEP 9910, 008 (1999). 

B. Zwiebach, JHEP , 0009, 028 (2000). 

J.A. Minahan and B. Zwiebach, JHEP 0103, :038 (2001). 
A. Sen, Int. J. Mod. Phys. A 20, 5513 (2005). 
M. C. Bento, O. Bertolami, and A. A. Sen, Phys. Rev. D 
67, 063511 (2003). 

A. T. Avelar, D. Bazeia, L. Losano, and R. Menezes, 
Euro. Phys. J. C 55, 133 (2008). 

R. Jackiw and C. Rebbi, Phys. Rev. D 13, 3398 (1976). 



